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ABSTRACT 

, We report the detection of CO (1-0) line emission in the bright, lensed star- forming 

galaxies SPT-S 233227-5358.5 (z = 2.73) and SPT-S 053816-5030.8 (z = 2.78), using 
the Australia Telescope Compact Array (ATCA). Both galaxies were discovered in 
a large-area millimeter survey with the South Pole Telescope (SPT) and found to 
be gravitationally lensed by intervening structures. The measured CO intensities im- 
ply galaxies with molecular gas masses of (3.2 ± 0.5) x 10 10 (^/15) _1 (Aco/0.8) M Q 
and (1.7 ± 0.3) x 10 10 (^/20)- 1 (A CO /0.8) M©, and gas depletion timescales of 4.9 x 
10 7 (Aco/0-8) yr and 2.6 x 10 7 (Aco/0-8) yr, respectively, where \i corresponds to the 
lens magnification and Aco is the CO luminosity to gas mass conversion factor. In the 
case of SPT-S 053816-5030.8, we also obtained significant detections of the rest-frame 
115.7 and 132.4 GHz radio continuum. Based on the radio to infrared spectral energy 
distribution and an assumed synchrotron spectral index, we find that 42 ± 10% and 
55 ± 13% of the flux at rest-frame 115.7 and 132.4 GHz arise from free-free emission. 
We find a radio-derived intrinsic star formation rate (SFR) of 470 ± 170 Mq yr -1 , con- 
sistent within the uncertainties with the infrared estimate. Based on the morphology 
of this object in the source plane, the derived gas mass and the possible flattening of 
the radio spectral index towards low frequencies, we argue that SPT-S 053816-5030.8 
exhibits properties compatible with a scaled-up local ultra-luminous infrared galaxy. 
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1 INTRODUCTION 

Submillimeter/millimeter blank field surveys have uncov- 
ered a significant population of dusty star-forming galax- 
ies (DSFGs) at cosmological distances. These galaxies are 
found to harbor intense star formation activity, with typical 
SF Rs > 500 Mm yr~ , and to be generally located at z > 
1 dBlain et al.l 120021: IChapman et ell 120051 ; lLagache et"al] 



l2005l ; IVieira et al.ll2013l ; IWeifi et al.ll2013f l. Their observa- 
tional properties suggest they are the likely progenitors of 
local massive early type galaxies, and they appear to trace 



large galaxy overdensiti e s at high redshifts dBrodwin et al 



2008 



Viero et all 



20ld ; ICapak et all 



200S| : iDaddi et al.l 120091; Aravena et al 



201ll ; lAmblard et al.ll201lh 



The interstellar medium (ISM) plays a critical role in 
shaping the overall evolution of galaxies as it constitutes the 
primary reservoir of material for star formation. The con- 
tinuous interaction between stars and the ISM governs the 
overall and chemical composition of galaxies. The molecular 
gas in the ISM is mostly composed of H2 and He, with a 
minor fraction (~ 1%) of other molecules and dust. How- 
ever, H2 is difficult to detect in the cold molecular ISM. The 
CO molecule, particularly in its lowest rotational transition 
J — 1 — 0, represents the best trac er of H2, spec ifically its 
mass and spatial distribution (e.g., IOmontll2007l ). 

Large efforts have been devoted to investigating the 
physical properties of DSFGs. However, due to their faint- 
ness and large cosmological distances, only a handful have 
been targeted in a sufficiently large number of molecular 
gas emission lines (e.g., CO) to allow for a full character- 
ization of their ISM. Follow-up CO studies have been fo- 
cused on the observation of J ^ 3 CO transitions. How- 
ever, CO (1-0) observations in a few unlensed DSFGs in- 
dicate that the gas masses are 2x larger and up to 3x 
more spatially extended than expected based on J 3 CO 
transitions, heavily altering our understanding of these sys- 
tems (e.g., Papadopoulos fc Ivisonl 20021; Greve et alj 2003; 



terns (e.g.. | rap adopoulos & ivison ZUUZ: l^revc et al J I AID J 
lHarris et al.ll2010l ; llvison et al.ll201ll ; iRiechers et alj|201lf > 
This highlights the need for low-J CO observations. 

Strongly lensed systems, including IRAS F10214+4722, 
the Cloverleaf, the Cosmic Eyelash and SMM 02399-0136, 
have proven to be a boon for studies of far-infrared 



and millimeter line studies at 



redshift 



Solomon fc Vanden Bout 2005: 


Brown & Vanden Bout 


1991 


; iRowan-Robinson et al. Il991; 


Brown & Vanden Bout 


1992 


; Barvainis et al. 


1992J; 


Ivison et al.ll 19981; iFraver et al. 


1998 


; Swinbank et al. 


2010 


). The recent discovery of a 



population of gravitationally lensed, highly magnified, 
star-forming galaxies in the > 100 square degree millimeter 
and far-I R surveys carried out b y the South Pole Telescope 
(SPT; ICarlstrom etafl l201ll) and the Herschel Space 
Observatory ( Pilbratt et al.l 120101 ) has provided a remark- 
able opportunity to study in great detail the process of 
star-formation and conditions of the IS M in this cosmolog- 
ically important po pulation of galaxies (|Vieira et alj|201ol ; 
iNegrello et afll201Ch . 

We have started a systematic follow-up campaign with 
the Australia Telescope Compact Array (ATCA) to study 
the molecular gas content in a sample of gravitationally 
lensed star-forming galaxies through the observation of 
J ^ 2 CO emission. These galaxies were found in the 1.4- 
3 mm wavelength sky maps produced by the SPT Sunyaev- 



Zel'do vich effect survey a nd have accurate spectroscopic red- 
shifts (|Greve et al.ll2012h . In this work, we report pilot ob- 
servations of the CO(1-0) line emission in two such sources. 
Hereafter, we adopt a ACDM cosmology with Ho = 71 km 
s" 1 Mpc" 1 , fl M = 0.27, and fi A = 0.73. 



2 OBSERVATIONS 
2.1 SPT sources 

The two targets of this work, SPT-S 233227-5358.5 
and SPT-S 053816-5030.8 (hereafter SPT 2332-53 and 
SPT 0538-50), wer e the first dusty-spectrum SPT sources 
l|Vieira et al.il201CJ ) to have spectroscopic redshift measure- 
ments. In both cases, a redshift determination was obtained 
through the identification of CO and [CI] line emission 
in wide-band submillimeter spectra obtained with the Z- 
spec instrument on the Atacama Pathfinder Experiment 
(APEX) 12m telescope, and confirmed with optical spec- 
troscopy obtained with the Very Large Telescope (VLT). 
The redshifts derived for the lensed objects are z = 2.738 
and z = 2.783 for SPT 2332-53 and SPT 0538-50, respec- 
tively l|Greve et al.ll2012l ). Both galaxies are associated with 
faint Spitzer Infrared Array Camera (IRAC) 4.5/xm sources 
(Fig. 1), and are detected with the APEX Large Bolometer 
Camera (LABOCA) and the Submillimeter APEX Bolome- 
ter Camera (SABOCA ) at 870^m and 350/im, respectively 
(see iGreve et all l2012h . 

As shown by the 850^im and 350/im images (Fig. [1] in 
IGreve et al1l2012h . the submillimeter emission in SPT 2332- 
53 is resolved into six major components, three of which 
are detected at 350pim (A, B and C; Fig. [lj. A massive 
foreground galaxy cluster at z = 0.403 is identified as the 
gravitational lens for this source l|Greve et alj|2012ft . A full 
analysis of the lensing structure of this object will be pre- 
sented in Vieira et al. (in preparation). That work leads to 
a magnification factor in the range fj, = 10 — 20, and thus 
we adopt n = 15 ± 5 throughout. 

SPT 0538-50 is unresolved in the LABOCA and 
SABOCA images. Recent ALMA 850 /im sub-arcsecond 
imaging of this source shows an Einstein ring around a fore- 
ground galaxy at z = 0.4 43, confirming the g ravitational 
lensed nature of this source l|Hezaveh et al.l2013li. A detailed 
lens m odel gives a magnification /i = 20 ± 4 ( Hezaveh et alj 
l2013h . Unless otherwise stated, the quoted uncertainties in 
the magnification factors for both sources are not propa- 
gate d to the obs er ved qu antities throughout. 

IGreve et all (|2012T) find observed IR luminosities 
Lir(8 - 1000/mi) = (6.5 ± 1.6) x 10 13 L Q for both sources. 



However, recent Herschel far-IR data leads to an improved 
observed IR luminosity of (8 .8 ± 1.3) x 10 13 Lq for SPT 
0538-50 (|Hezaveh et al.l 120131 ). The revised value for the 
IR luminosity of SPT 2332-53 is w e ll wit hin uncertainties 
with the estimate from lGreve et all l|2012h . and we use this 
value throughout. From this, we adopt intrinsic IR lumi- 
nosities of (4.3 ±1.1) x 10 12 (At/15)" 1 Lq and (4.4 ±0.65) x 
10 12 (/i/20) _1 L Q for SPT 2332-53 and SPT 0538-50, respec- 
tively. The estimated total IR luminosities are derived from 
fitting a single-temperature grey-body model to the data at 
Arcst > 50 jj,m, and t hus do not accoun t for a possible hot 
dust component (see IGreve et alj |2012. for details). While 
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Figure 1. CO(l-O) and 35 GHz continuum detections toward SPT 2332-53 and SPT 0538-50. The Spitzer IRAC 4.5/im image is shown 
in greyscale with blue contours overlaid representing the CO(1-0) line emission (dirty image) averaged over 350 km s -1 and 220 km s -1 , 
respectively. These contours are in steps of ±ler, starting at ±2cr, with a = 0.32 mjy beam -1 and a = 0.44 mjy beam -1 , respectively. 
Red contours represent the 35 GHz continuum emission (dirty image), and are shown in steps of +lcr, starting at +2c, with <r = 31 fijy 
beam -1 and a = 21 fijy b eam -1 , respectively . For SPT 2332-53, the location of the SABOCA sources are shown as red plus symbols, 
following the notation from lGreve et al] J2012h . The axes show the offset, in arcseconds, with respect to the phase tracking center. 



such hot dust component would increase the actual IR lumi- 
nosities of our targets, such emission would likely arise from 
AGN and thus would not strictly trace the star formation 
in the host galaxies, being irrelevant for the analysis below. 



2.2 ATCA CO data 

We used ATCA in the H75 compact-array configuration 
to observe the redshifted CO(l-O) emission line (v tes t = 
115.271 GHz) in the two star-forming galaxies SPT 2332-53 
and SPT 0538-50. The observations were performed between 
UT 2012 July 27-30 under good weather conditions with 5 
working antennas, and the sources were observed for ~10h 
each. 

We used the Compact Array Broadband Backend 
(CABB) configured in the lM-0.5k mode, which leads to 
a bandwidth of 2 GHz per correlator window with 1 MHz 
per channel of spectral resolution. In both cases, one of 
the windows was tuned to 35000 MHz to obtain a mea- 
surement of the continuum emission. For SPT 2332-53 the 
other window was tuned to 31412 MHz, and for SPT 0538- 
50 this was tuned to 30600 MHz. The phase tracking cen- 
ters were a 2 ooo = 23 h 32 m 27.8 3 , S 20 oo = -53°58'37", and 
a 2 ooo = 05 h 38 m 16.8 s , 8 20 oo = -50°30'52", respectively, 
based on the prior positions obtained with LABOCA at 
870^im. 

The bright quasars 2355-534 (S35GHZ ~ 1-8 Jy) and 
0537-441 (535GHz ~ 6.2 Jy) were observed every 6-7 min for 
gain calibration, and every ~ lh for pointing calibration. 
The sources 1921-293 (~ 14 Jy) and 0537-441 were used 
for bandpass calibration, and Uranus and 1934-638 were 
used as amplitude calibrators. The software package Miriad 
ijSault et al.l Il995l ) an d the Common Astro nomy Software 
Applications (CAS A; iMcMullin et all 120071 ) were used for 



editing, calibration, and imaging. Absolute flux calibration 
is estimated to be accurate to within 15%. 

Flagging of visibilities was performed for the case of 
SPT 0538-50, where some channels fell close to the edge of 
the Q-band short frequency limit at 30 GHz. The visibilities 
were inverted using natural weighting yielding synthesized 
beam sizes for SPT 2332-53 and SPT 0538-50 of 21.8" x 16.8" 
and 22.2" x 13.9" (at ~ 31 GHz), and 19.2" x 14.8" and 
18.0" x 14.2" (at 35 GHz), respectively. The final data results 
in an rms of 0.70 mjy beam -1 per 4 MHz (39 km s -1 ) chan- 
nel at ~ 31.4 GHz and 0.60 mjy beam -1 per 2 MHz (19.7 
km s -1 ) channel at ~ 30.6 GHz for SPT 2332-53 and SPT 
0538-50. We reach continuum rms levels of 40 /iJy beam - 
and 35 /iJy beam -1 at 31.4 and 35 GHz, respectively, for 
SPT 2332-53, and 35 /xJy beam -1 and 20 /iJy beam -1 at 
30.6 and 35 GHz, respectively, for SPT 0538-50. 



2.3 ATCA low frequency data 

ATCA observations at 2.2, 5.5, and 9.0 GHz of SPT 2332-53 
and SPT 0538-50 were made between UT 2012 January 23- 
27 using the CABB in the lM-0.5k mode (see above). The 
observations were performed in the most extended ATCA 
configuration, 6A, with 6 working antennas. The data were 
edited, calibrated and imaged using the Miriad package. 
Data affected by known radio frequency interference (RFI) 
or with bad visibilility ranges were flagged accordingly. We 
estimate an absolute calibration uncertainty of ~ 5% at 2.2 
and 5.5 GHz, and 10% at 9.0 GHz. We inverted the visi- 
bilities using natural weighting leading to beam sizes of 
~ 7.7" x 5.2", 3.2" x 2.1" and 2.0" x 1.3" at 2.2, 5.5 and 
9.0 GHz, respectively. For SPT 2332-53, we reach rms noise 
values of 200, 500, and 100 ^iJy beam -1 at 2.2, 5.5, and 
9.0 GHz, respectively. For SPT 0538-50, we reach rms noise 
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values of 200, 100, and 75 fiJy beam 1 at each frequency, 
respectively. 



3 RESULTS 
3.1 SPT 2332-53 

The CO (1-0) line toward SPT 2332-53 is well detected and 
spatially resolved into two main peaks separated by ~ 30" 
(Fig. [TJ. The CO profiles of both peaks (east and west) 
are coincident in frequency and line shape suggesting that 
they arise from the same lensed background ob ject (Fi g. [3). 
A com parison with the SABOCA image from|c reve et al.l 
(2012) indicates that the CO distribution corresponds to the 
blended IR components A, B, and C, following the nomen- 
clature adopted by Greve et al. 

The continuum emission was measured from the 35 GHz 
correlator window and the line-free channels at ~31.4 GHz. 
We find significant detections at both frequencies from a 
source located at ~ 15" to the north of the CO emis- 
sion, with flux densities S31.4GHZ = 310 ± 40 fiJy and 
S35GHZ — 300 ±35 /iJy (Fig. [TJ This source is identified with 
the brightest cluster galaxy (BCG) in the foreground struc- 
ture. It is detected at 3.4 and 4.6 /im with Spitzer IRAC, 
though not at 12 and 22 fim with the Wide-field Infrared 
Survey Explorer (WISE). This source is also detected in the 
5.5 GHz and 9.0 GHz ATCA images. The non-detection of 
continuum at the CO position yields a limit to the inte- 
grated ~ 35 GHz flux density of the background source of 
< 220 /iJy (3a with respect to the zero level). 

From Gaussian fitting to the CO profile of the combined 
emission from both components (Fig. [2]), we find the fluxes 
and line-widths listed in Table Q] The total line intensity 
is obtained by creating a moment-0 image, integrated over 
~ 600 km s~ x around the line (Table[TJ). 



4 ANALYSIS 



4.1 The nature of the SPT sources 



The results of the source reconstruction of SPT 2332-53 do 
not allow us to discern whether this object corresponds to a 
disk galaxy or merger (Vieira et al in prep.). The CO profile 
shape, however, does not show evidence for a double-peak 
structure (Fig. [2J as would be exp ected for a disk geo metry 
extended on scales of ~ 5 — 10 kpc (Da ddi et al.|[20 10) or for 
an early-stage merging system (e.g., Greve et al.l2005h . This 
suggests the molecular gas is distributed over a relatively 
compact disk or a compact spherical distribution, as in an 
advanced merger stage, with most of the gas within a radius 
~ 1 - 2 kpc. 

In the case of SPT 0538-50, the lens reconstruction 
shows that this sour ce is composed of two components (see 
iHezaveh et all 120131 ): a compact source with a radius of 
0.5 ± 0.1 kpc, and an extended source with a radius of 
1.6 ±0.3 kpc. The two components appear to spatially over- 
lap, leaving it unclear whether the system corresponds to 
two interacting galaxies, or a clumpy disk with a bright com- 
pact source. However, the fraction of the ALMA 870 fim flux 
distributed between both components, ~ 70% of the total 
Lir for the compact s ource and ~ 30% for the extended one 
(|Hezaveh et al.ll2013h . is in good agreement with the CO 
luminosity budget. The brighter CO component (SPT 0538- 
50blue) carries a luminosity L' co = 1.49 x 10 10 (fi/20)- 1 K 
km s _1 pc 2 or ~70% of the total CO luminosity, while the 
fainter CO component (SPT 0538-50red) has a CO luminos- 
ity L' co = 6.1 x 10 9 (m/20)" 1 K km s" 1 pc 2 or ~30%. This 
suggests that the compact reconstructed source could be 
identified with the bright CO line component, while the ex- 
tended source could be associated with the fainter, broader 
CO line component (Fig. [2J. However, this IR to CO line 
association remains uncertain. The only way to confirm this 
is through spatially resolved CO spectroscopy. 



3.2 SPT 0538-50 

The CO (1-0) line toward SPT 0538-50 is also significantly 
detected, however at the spatial resolution of these observa- 
tions, the source is unresolved (Fig. [T]). The spectrum sug- 
gests the existence of two CO line peaks (Fig. [2|. Fitting 
a two component Gaussian profile to the spectrum, we find 
the bright (SPT 0538-50blue) and weak (SPT 0538-50red) 
CO components separated by 350 ± 50 km s _1 in velocity. 
The fainter CO component carries «30% of the CO inten- 
sity. Table[T]summarizes the measurements for these sources. 
We measure the total CO line intensity from the moment-0 
image, integrated over 710 km s" 1 containing both CO line 
peaks (SPT 0538-50; Tabled]). 

Continuum emission at 30.6 and 35 GHz is also detected 
(Fig. [TJ. In this case, the continuum peak position coincides 
within l'.'O of the location of the CO peak. The origin of this 
emission is discussed in the next sections. Continuum emis- 
sion is also detected at 2.2 and 5.5 GHz, with flux densities 
of 52.2GHz = 650 ±200 /iJy and Ss.sghz = 550±100AtJy, re- 
spectively. This source is not detected in our 9.0 GHz data, 
resulting in a 3<r upper limit of 225 /iJy (respect to the zero 
level) . 



4.2 Gas and dynamical masses 

The molecular gas mass, Af gas is commonly computed from 
the CO (1-0) line luminosity using the relation M gas = 
A'co-LcO! where Xco is the CO luminosity to gas mass 
conversion factor. Its value c an vary from galaxy to galaxy, 
and even within galaxies (e.g., G enzel et al.ll2012r ). We adopt 
Xco = 0.8 Mq (K km s _1 pc 2 ) -1 , a s found for local ul- 
traluminous infrared galaxies (ULIRGs; lDownes fc Solomon! 
1998). While this value is typically used for dusty starburst 
galaxies at high redshift, higher values have b een recently 
found for distant disk galaxies (Xco ~ 3.6; Dad di et al.l 
l201Ct ). Using a larger value for Xco would produce gas 
masses 4.5 times larger than quoted here. From the CO lu- 
minosities found in the previous section, we find M gas = 
(3.2 ±0.5) x 10 10 (^/15)- 1 (A CO /0.8) M Q and M gas = (1.7± 
0.3) x 10 10 (/i/20)- 1 (A C o/0.8) M© for SPT 2332-53 and 
SPT 0538-50, respectively. 

We compute dynamical masses for our sources based 
on the CO line profiles. For a disk gas distribution, the 
dynamical mass of the system is given by Md yn ~ 2.3 x 
10 5 Avp WHM R sin _2 (i) Mq , where Avfwhm is the ob- 
served CO line FWHM in units of km s _1 , R is the disk 
radius in kpc, and i is the inclination angle. For a merger 
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Table 1. Observed CO properties 



Source 


zco a 


Sco b 


"FWHM 6 


^CO 1-0 c 


Tl d 


S31GHZ e 


S35GHZ e 






(mjy) 


km s — 1 


(Jy km s — 1 ) 


(xlO 10 in/no)' 1 K km s" 1 pc 2 ) 


(MJy) 


(MJy) 


SPT 2332-53t 


2.7256(2) 


4.66 ± 0.51 


342 ± 42 


1.70 ±0.25 


3.94 ±0.58 


< 120 


< 120 


SPT 0538-50blue* 


2.7854(1) 


3.81 ± 0.36 


210 ±34 


0.85 ±0.08 


1.49 ±0.14 






SPT 0538-50red* 


2.7897(6) 


1.18 ±0.27 


281 ± 80 


0.35 ±0.18 


0.61 ±0.31 






SPT 0538-50 


2.7855(1) 






1.20 ±0.20 


2.16 ±0.36 


171 ±35 


130 ± 20 



t Integrated values over all CO components are listed. 

t Here, SPT 0538-50blue and SPT 0538-50red are the brighter and fainter CO peaks, respectively. 

a CO redshift. For SPT 0538-50, the redshift is computed as the weighted average between the redshifts of SPT 0538-50blue and red. 
b CO line peak flux and line full-width half maximum (FWHM) velocity. 
c Integrated CO line intensity (Ico = f Scodv). 

d CO line luminosity. The magnification factor fio is assumed to be 15 and 20 for SPT 2332-53 and SPT 0538-50, respectively. 
e Continuum fluxes at ~ 31 GHz and 35 GHz. 
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Figure 2. ATCA spectra of the CO(1-0) line emission toward SPT 2332-53 and SPT 0538-50. Thei r spectra are shown at 39 km s" 1 (4 
MHz) and 20 km s — 1 (2 MHz) resolution, respectively, and centered at the redshift determined bv lGreve et al. I (I2012T) . A solid red line 
shows a two-component Gaussian fitting to the CO profiles in each case. The inset plot for SPT 2232-52 shows the comparison of the 
CO profiles for the eastern and western components. 



model, Awfwhm represents half the physical separation in 
velocity of the component CO lines, and R is half the separ a- 
tion between components (|Solomon &: Vanden Boutl [2005). 



Given the poor constrains on the geometries of both 
sources, we derive dynamical masses by adopting disk ge- 
ometries. For SPT 2332-53, we use R = 2 kpc and an aver- 
age inclination ang le i ~ 57° (|Law et al.ll2009r i, which yields 
M dyn = 7.8 x 10 10 Mq. This implies a gas fraction of ~ 0.4 
for jj, — 15 and the adopted Xco value. For SPT 0538-50, we 
use R — 1.6 kpc (from the lens model), a CO profile FWHM 
of ~ 500 km s" 1 (from the full CO line) and an average in- 
clination i ~ 57°, yielding M dyn = 1.3 x 10 11 M©. This leads 
to a gas fraction of ~0.1 for /1 = 20 and the adopted Xrjo- 
We remark that since the actual source geometries are un- 
known, the estimated dynamical masses and gas fractions 
are uncertain and should be taken with caution as they are 
given for reference only. 



4.3 Star formation efficiencies and depletion 
timescales 



The star formation efficiency (SFE) can be defined as SFE 
= Lir,/L co in units of Lq (K km s _1 pc 2 ) -1 . From this, 
we obtain SFEs of 110 and 205 L e (K km s" 1 pc 2 )' 1 for 
SPT 2332-53 and SPT 0538-50, respectively. The gas deple- 
tion timescale can be defined as ta ep ~ M gas / SFR. We adopt 
SFR (M© yr" 1 ) = 1.5 x IP" 10 L m (L^) jMurphv et all 
l201ll) for consistency with iGreve et all l|2012l ). This yields 
gas depletion timescales of ~ 4.9xl0 7 (/x/15) _1 (Xco/0.8) vr 
and ~ 2.6 x 10 7 (^/20) _1 (Xco/0.8) yr for SPT 2332-53 and 
SPT 0538-50, respectively. While values for the depletion 
timescal es for DSFGs and l ocal ULIRGs are typically < 10 8 
yr (e.g.; IGreve et al.ll2005l ; [Solomon fc Vanden BoutJ|2005l ). 
distant "normal" star forming disk galaxies w ere found to 
have depletion timescales o f ~ (0 .3—2) x 10 9 yr (Dad di~et al.l 
l20ld ; iTacconi et ai1l201(j 120131 '). The depletion timescales 
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Figure 3. The ALMA 345 GHz high -resolution map for 
SPT 0538-50 is shown in the background l lHezaveh et alj|2013t) 
with red contours overlaid representing the ATCA 5.5 GHz radio 
map. Contour levels are in steps of ±lc7, starting at ±2cr with 
<t — 35 (j,Jy beam -1 . The spatial coincidence between the sub- 
millimeter and radio emission confirms that the latter arises from 
the background lensed object. 

for our sources appear to be consistent with typical DSFGs 
and local ULIRGs. 



4.4 SPT 0538-50: Free-free, synchrotron, and dust 
emission 

SPT 0538-50 is additionally detected at 2.2 and 5.5 GHz by 
ATCA. At z = 2.783 the 2.2, 5.5, 30.6 and 35 GHz obser- 
vations sample 8.3, 20.8, 115.7 and 132.4 GHz in the rest- 
frame, respectively, providing a long lever-arm to measure 
the radio spectral index of the source. 

Fig.[3]shows the recently obtained ATCA 5.5 GHz con- 
tinuum image toward SPT 0538-50. The image structure 
matches the position a nd structure of the Ein stein ring found 
in the ALMA images l|Hezaveh et al.lfeoist ). clear evidence 
that the bulk of the radio emission comes from the lensed 
background source. We thus rule out the foreground galaxy 
contribution to the radio emission based on the morphology 
shown by ALMA and ATCA. 

In the following section, we analyze the radio spectral 
properties of SPT 0538-50, and adopt the rest-frame fre- 
quencies for clarity. We define the flux density S v in terms 
of the observed frequency v as S v oc v a , with a the radio 
spectral index. The data seem to suggest that the radio spec- 
trum may steepen between rest-frame 8.3 and 132.4 GHz as 
the observed spectral index between 8.3 and 20.8 GHz is 
-0.18 ±0.39, and -0.76 ±0.11 between 20.8 and 132.4 GHz 
(see Fig. [4}. In fact, subtracting off an estimate of the ther- 
mal dust contribution to the 132.4 GHz flux density suggests 
an even steeper intrinsic (free-free ± synchrotron) radio 
spectral index between 20.8 and 132.4 GHz of -0.89 ± 0.14. 
Given the large uncertainties on the photometry, we do not 
include an additional free parameter for spectral curvature 
in our fit to the radio spectrum, and simply assume a sin- 
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Figure 4. The radio-to-IR spectrum of SPT 0538-50 fit by a 
combination of modeled radio and IR spectra (solid line). Pho- 
tometry is shown with error bars, along with an upper limit at 
frest =25 GHz. The IR dust emission is given by a grey-body 
model (dotted line). The radio spectrum is modeled by 2 compo- 
nents, non-thermal synchrotron emission (long-dashed lines) and 
thermal, free-free, emission (dot dashed line) 

gle power law without a break. However, we consider and 
discuss the physical implications of the observed spectral 
flattening at lower frequencies below in Section [4.4.31 

4-4-1 Fitting the radio-to-IR spectrum: free-free emission 
and an independent SFR estimate 

To estimate the amount of free-free emission at rest-frame 
115.7 and 132.4 GHz, which can then be used to calculate a 
SFR, we fit the full radio-t o-IR spectrum. For simplicity, we 
use the grey-body fit from lGreve et al.l (|2012h to model the 
thermal dust component, which has a fixed emissivity spec- 
tral index f} = 2. To fit the radio spectrum, we vary a combi- 
nation of free-free (thermal) and synchrotron (non-thermal) 
radio emission components which scale as oc c° T and 
S^ T oc u aNT , respectively, where ckt = —0.1 and <»nt are 
the free-free and synchrotron spectral indices, respectively. 

Before varying the free-free and synchrotron compo- 
nents to fit the observed radio spectrum, we first need to 
assume a non-thermal spectral index. We assume a syn- 
chrotron spectral index of cint = —1.0, which is consis- 
tent with the steepest part of the radio spectrum measured 
between 20.8 and 132.4GHz (see above). This value is in 
agreement with an estimate for the non-thermal spectrum 
based on the calculation of energy losses to cosmic ray elec- 
trons in this system arising from synchrotron radiation, in- 
verse Compton scattering, ionization, brems strahlung, and 
escape through an empirical prescription (see lMurphv||2009l . 
for details). Furthermore, this value is consistent with the 
synchrotron spectral index of —1.05 ± 0.0 5 for two dusty 
star-fo rming galaxies at z ~ 2.5 found by I Thomson et al.l 
|2012l ). Setting qnt = —1.0, our best fit model to the radio- 
to-IR data suggests a free-free radio fraction at 115.7 and 
132.4 GHz of 42 ± 10% and 55 ± 13%, respectively (see Fig. 
©. 

Taki ng these value s and the free- free SFR calibration 
given in iMurphv et all (|2012T l results in a SFR = 470 ± 
17OM0yr _1 . This is consistent with the IR luminosity- 
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derived estimate of SFR = 650 ± 100 M© yr" 1 , which is 
again calculated using the calibration given in Murphv et al. 
[|201lL calibrated for a common IMF (i.e.. lKroupal | 200lh . 
and consistent with the estimate from lGreve et alTl|2012l i. 
We note that such calibrations between IR luminosity and 
SFR are not lik e ly accurate to bet ter than a factor or 2 
l|Kennicuttlll998l ; lOmont et ai1l200ll ). and this uncertainty 
is not included in the error above on the IR-derived SFR. 



4-4-2 Physical considerations in the fitting 

While we have assumed a physically motivated value for the 
non-thermal spectral index, it is worth noting that this as- 
sumed value is marginally inconsistent with a naive measure- 
ment of a from our data. Using an ordinary-least-squares fit 
to the four radio detections, weighted by the uncertainties 
of the flux density measurements, results in a radio spec- 
tral index (before correcting for thermal dust emission) of 
ct = —0.65 ± 0.10. Assuming this value for the non-thermal 
radio spectral index, the fit results in a 132.4 GHz free-free 
fraction of ~2%. Similarly, using a dust emissivity index 
j3 = 2 and a radio spectral index ~ 0.75 l|lbar et al.ll201Cf ) 
would result in a free- free fraction of ~20%. 

For any physical configuration in which most of the IR 
emission is powered by star formation, such a small frac- 
tion of free-free radio emission is highly unlikely. A situa- 
tion without free-free emission could arise for cases where: 
(a) the emission is purely powered by an AGN; (b) there is 
a suppression of Lyman continuum photons, and thus free- 
free emission, through an upper mass cutoff of the IMF at 
~ 20 Mq while still providing synchrotron emission from su- 
pernovae (~ 8 — 20 Mq); and (c) the absorption of ionizing 
photons by dust in the actual starburst. However, the galaxy 
shows a rather normal, albeit slightly higher, IR/radio ra- 
tio (see Section f4.4.3|) . which is difficult to keep constant if 
there are variations in the IMF or in the presence of a pow- 
erful AGN. In the case where dust successfully competes 
with neutral hydrogen for ionizing photons, we still expect 
a significant free- free component to be present and compati- 
ble with the amou nt of ongoing star for mation inferred from 
the IR spectrum (|Murphv et alj|2012h . A similar situation 
has been observed in local (U)LIRGs, where the synchrotron 
spectrum is fo und to be increasing l y depressed at hig h fre- 
quencies (e.g.. IClemens et al1l2008l ; iLerov et alj|2012h . per- 
haps arising from a modified electron injection spectrum in 
such dense starbursts. 



4-4-3 Spectral Flattening Due to a Deeply Embedded 
Compact Starburst 

The radio observations suggest that the radio spectrum 
flattens towards lower frequencies. This flattening may be 
the result of the radio continuum emission becoming opti- 
cally thick at lower frequencies, a s observed in local com- 
pact ULIRGs (jCondon et al.l fl991V Indeed, the lens model 
of SPT 0538-50 indicates the presence of an IR-dominant, 
compact source. Such a spectral flattening at lower ra- 
dio frequencies typically causes sources to have IR/radio 
(1.4 GHz) ratios that are la rger than the aver age value 
among star- forming galaxies |Condon et al.|[l99 lh . Extrap- 
olating a 1.4 GHz flux density using the observed rest-frame 



8.3 GHz flux density along with the spectral index mea- 
sured between 8.3 and 20.8 GHz suggests that the galaxy 
has a logarithmic IR/1.4GHz ratio qm, — 2.843- This im- 
plies a linear ratio nearly 60% (l<r) larger than the average 
value of qiR = 2.64 ±0.02 found for local star-forming galax- 
ies with 1.4 GHz luminosities in the range ~ 10 19 - io 24 - 5 
W Hz- 1 (|Bellll20"03h . 

Assuming that the spectral flattening is in fact the re- 
sult of the source being powered by an optically-thick star- 
burst, we can us e local relation s to in fer properties of the 
galaxy. Recently iMurphv et ail (|2013h found a correlation 
between the radio spectral index of z ~ (U)LIRGs and 
their distance from the star-forming main sequence (i.e., 
specific SFR defined as sSFR = SFR/M*). That is, galaxies 
having flatter spectral indices, most likely due to harboring 
compact starbursts, which are optically thick at low radio 
frequencies, have increasingly higher sSFRs. Based on this 
correlation, and assuming an evolution in the sSFR of galax- 
ies from 2 ~ to z ~ 3, where the typical sSFR at z ~ 3 
is 4.5 Gyr -1 l|Karim et alj l201ll ). the spectral index mea- 
sured between 8.3 and 20.8 GHz suggests that the sSFR of 
SPT 0538-50 is ~ 10.2 ± 2.9 Gyr -1 , more than a factor of 
~2 above the z ~ 3 main sequence for star-forming galaxies 
l|Karim et al.ll201lT ). Taking the above free- free estimate for 
the SFR suggests this system has M* = (4.6±2.1) x 10 10 M . 



5 DISCUSSION AND SUMMARY 

We have obtained sensitive CO(1-0) and radio continuum 
ATCA observations toward two bright DSFGs discovered 
in millimeter survey data taken with the SPT. This cor- 
responds to a pilot study of a larger sample of millimeter 
bright SPT sources with available spectroscopic redshifts. 

Both systems are found to be rich in molecular gas 
with masses ~ 2 x 10 10 M Q , gas fractions of - 0.1 - 0.4 
and SFEs of ~ 100 - 200 L© (K km s _1 pc 2 ) -1 , for the 
adopted values of the magnification (fj,) and CO luminosity 
to gas mass conversion factor (Ago). Taken together, the 
observed CO properties are compatible with those found 
in other lensed (and unlensed) submillimeter-se lected galax- 
ies based on previous CO (1-0) measureme nts (|Fraver et al.l 
|20U| ; llvison et al.ll201ll ; lHarris et al]|2010h . 

In both cases, we find ~ 31 and 35 GHz continuum de- 
tections in the field (rest-frame 115.7 and 132.4 GHz). While 
in the case of SPT 2332-53 the continuum emission comes 
from the foreground BCG galaxy, in the case of SPT 0538-50 
the radio emission is clearly identified with the CO source, 
as shown in Fig. [3] Assuming a non-thermal spectral index 
of qnt = —1.0, our modeling of the radio to IR spectrum 
of the SPT 0538-50 suggests that 42 ± 10% and 55 ± 13% 
of the emission at rest-frame 115.7 and 132.4 GHz, arise 
from a free- free component. This allows us to derive a free- 
free SFR value of 470 ± 170 M© yr -1 , which is consistent 
within the unce r taintie s with the IR-derived SFR. Recently, 
iThomson et al.1 (|2012h detected CO (1-0) line emission m 

1 The latter is defined as the ratio of the total IR luminosity 
(8 — lOOO^tm) to the radio power: 

, ( £lR \ . / 5i.4GHz \ 

qm = l0g10 U75X1012W m-' ) " 10810 Um-' Hz"* J ' 
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three submillimeter-selected galaxies at z ~ 2.5 — 2.9, and 
marginally detected the rest-frame 115 GHz continuum from 
two of these galaxies. Their modeling of the radio spectra 
indicates free- free contri butions to the rest-fra me 115 GHz 
in the range - 30 - 40% l|Thomson et al.ll2012T l. This is sim- 
ilar to our results, which are based on significant detections 
at rest-frame 115.7 and 132.4 GHz on SPT 0538-50. 

Analysis of the radio spectral curvature for this source 
suggests that the star formation activity is powered by a 
compact starburst; the spectral flattening towards lower 
radio frequencies is similar to what is observed for lo- 
cal (U)LIRGs harboring compact starbursts that become 
optically-thick at low radio frequencies. From this analy- 
sis, we estimate a stellar mass of (4.6 ± 2.1) x 10 10 Mq. 
It is interesting to note that almost all local ULIRGs are 
powered by merger-driven star formation and lie above the 
star-forming main sequence, which also seems to be the case 
for SPT 0538-50 with a sSFR of 10.2 Gyr" 1 compared to 
the main-sequence at z ~ 3 with an average sSFR of 4.5 
Gyr -1 . Furthermore, the ULIRG-like nature of this source 
is supported by the results from the lens model which in- 
dicate that SPT 0538-50 is composed of a compact IR- 
dominant component, which c omprises 70% of the IR lu- 
minosity (|Hezaveh et al. 20131). This is simila r to what is 
seen in local ULIRGs (|Brvant fe ScoviUdll999h . 

These results exemplify the use of deep low-J CO spec- 
troscopy in high-redshift star forming galaxies taking ad- 
vantage of the large bandwidths available with ATCA and 
the Karl Jansky Very Large Array (VLA). The advent of 
important samples of bright, gravitationally lensed galaxies 
from the SPT and Herschel surveys, make possible detailed 
studies of the ISM in DSFGs. While observation of the CO 
line emission enables a direct measurement of the gas and 
dynamical masses, modeling of the radio spectrum enables 
independent measurements of the SFRs and stellar masses 
without the need of near-to-far IR imaging, which generally 
requires space-based observatories. 
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